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ABSTRACT 

In this paper, we present a comprehensive analysis of star-forming galaxies (SFGs) at intermedi¬ 
ate redshifts (z^l). We combine the ultra-deep optical spectro-photometric data from the SHARDS 
survey with deep UV-to-FIR observations in the GOODS-N field. Fxploiting two of the 25 SHARDS 
medium-band filters, F687W17 and F823W17, we select [OH] emission line galaxies (FLGs) at 
2^0.84 and z^l.23 and characterize their physical properties. Their rest-frame equivalent widths 
(FWrf([OH])), line fluxes, luminosities, star formation rates (SFRs) and dust attenuation properties 
are investigated. The evolution of the BWrf([OH]) closely follows the SFR density evolution of the 
Universe, with a FWrf([OH])cx(l-|- 2 :)^ trend up to redshift z~l, followed by a possible flattening. 

The SF properties of the galaxies selected on the basis of their [OH] emission are compared with 
complementary samples of SFGs selected by their MIR and FIR emission, and also with a general 
mass-selected sample of galaxies at the same redshifts. We observationally demonstrate that the UVJ 
diagram (or, similarly, a cut in the specific SFR) is only partially able to distinguish the quiescent 
galaxies from the SFGs. The SFR-M* relation is investigated for the different samples, finding a 
logarithmic slope ~1, in good agreement with previous results. The dust attenuations derived from 
different SFR indicators (UV(1600), UV(2800), [OH], IR) are compared, finding clear trends with 
respect to both the stellar mass and total SFR, with more massive and highly star-forming galaxies 
being affected by stronger dust attenuation. 

Subject headings: galaxies: emission lines — galaxies: star formation — galaxies: high-redshift — 
galaxies: evolution — galaxies: general — galaxies: photometry 


1. INTRODUCTION 

The study of emission line galaxies (FLGs) has long 
been recognized as one of the most powerful and direct 
tools to investigate star-forming galaxies (SFGs) at dif¬ 
ferent redshifts. For local SFGs, the most common and 
reliable emission line used to investigate star formation 
(SF) is Ha (see e.g. Gallego et al. 1995; Tresse & Maddox 
1998; Perez-Gonzalez et al. 2003; Dominguez Sanchez et 
al. 2012). At increasing redshifts other lines and methods 
must be used because the Ha line runs out of the optical 
and near infrared (NIR) range, becoming more difficult 
to be observed. Present day NIR facilities allow to extend 
these local studies to higher redshifts using spectroscopy 
(see e.g. Tresse et al. 2002; Twite et al. 2012) or narrow- 
band photometry (see, e.g. Villar et al. 2008, 2011; Sobral 


et al. 2009b, 2011, 2013; Ly et al. 2011), although it re¬ 
mains hard and expensive, in terms of observing time, to 
obtain deep and wide surveys of line emitters with these 
instruments. 

One of the most popular SF indicators used at interme¬ 
diate/high redshift is the [OH] emission line at 3727 A. 
This line lies at optical wavelengths up to z^^l.b, so 
it has been widely use to detect and study SFGs for 
decades. The physical mechanism responsible for the 
production of this intense emission line (actually a dou¬ 
blet at 3726 A and 3729 A) has been investigated and well 
understood for a long time (Fddington 1927), but the 
relation between environmental conditions (density, tem¬ 
perature, metallicity, ionization state of the gas, ...) and 
the intrinsic emitted flux is still far from being completely 
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characterized. Nonetheless, [Oil] has demonstrated to be 
a very useful SFR indicator and empirical calibrations 
are commonly used to transform [Oil] luminosities into 
SFR (Gallagher et al. 1989; Kennicutt 1998; Hogg et al. 
1998; Jansen et al. 2001; Kewley et al. 2004; Yan et al. 
2006; Moustakas & Kennicutt 2006; Garn & Best 2010; 
Gilbank et al. 2010, 2011; Hayashi et al. 2013). 

The Survey for High-z Absorption Red and Dead 
Sources (SHARDS^, Perez-Gonzalez et al. 2013) is an 
ESO/GTG Large Program carried out with the OSIRIS 
instrument on the 10.4m Gran Telescopio Canarias 
(GTG). Our survey has obtained imaging data in the 
GOODS-N field through 25 medium-band filters cover¬ 
ing the wavelength range between 500 and 950 nm in 
sub-arcsec seeing conditions. The typical width of our 
filters is 15 nm for wavelengths bluer than 880 nm, and 
25-35 nm for the three reddest filters. In each filter, 
SHARDS is able to detect 26.5 mag sources at the do- 
level (at least). SHARDS allows us to study any iso¬ 
lated emission line falling within its wavelength coverage. 
More specifically, for the [OH] line, we can cover the in¬ 
terval from z^O.3 up to Z'^1.6. In this work, we focus 
on the analysis of ELGs selected using two of the twenty 
five SHARDS filters, the F687W17 and F823W17, cor¬ 
responding to ELGs at z~0.84 and z~1.23 respectively. 
We have chosen these two filters because they are repre¬ 
sentative enough of the whole dataset and they provide 
two samples of galaxies that are sufficiently separated in 
redshift (~1.5 Gyr in time) to allow an investigation of 
the possible evolutionary effects on the derived physical 
properties. A more extensive analysis including all the 
SHARDS filters will be presented in a future work. 

One of the aims of this work is to demonstrate the 
the power of the SHARDS medium-band data to select 
ELGs down to very faint continuum magnitudes and line 
fluxes (similar or even fainter that those characteristic of 
the deepest spectroscopic and narrow-band surveys) and 
study their physical properties in detail. Even though 
the equivalent spectral resolution of SHARDS (R^50) 
is smaller than the one reached by typical spectroscopic 
(e.g. TKRS, VVDS, DEEPS, HETDEX or MOSDEF 
surveys; Wirth et al. 2004; Le Eevre et al. 2005; Cooper 
et al. 2011; Adams et al. 2011; Kriek et al. 2015) or nar¬ 
row band surveys (Villar et al. 2008, 2011; Sobral et al. 
2009b, 2010, 2011, 2012, 2013), the depth and quality of 
the SHARDS data compensate this a priori disadvantage 
allowing us to study stellar mass complete samples down 
to M* ~10® Mgun at z^l (Hernan-Caballero et al. 2014). 

In fact, our medium-band survey can reach deeper 
magnitudes (down to R'^26-27) with respect to typi¬ 
cal spectroscopic survey (limited to R~24-25) by in¬ 
vesting substantially lower amount of observing time 
and overcoming the issues arising when dealing with 
slit apertures and the limited multiplexing of spectro¬ 
graphs. On the other hand, integral field unit (IFU) spec¬ 
trographs can also observe all the sources within their 
field of view (FOV), but this is normally very limited, 
while the SHARDS filters apply to the relatively wide 
OSIRIS@GTC field of view (i.e. ~ 7' x 8'). Further¬ 
more, ultra-deep medium-band selected SF galaxy sam¬ 
ples should be less affected by Malmquist bias with re- 

^ The SHARDS web page is available at: 
http: //guaix.fis.ucm.es/'^pgperez/SHARDS/ 


spect to spectroscopically observed samples which tend 
to favour highly SFGs. We also remark that system¬ 
atic effects are expected to be introduced in the study of 
emission-line luminosity functions by the typical selec¬ 
tion of targets in spectroscopic surveys, which are typ¬ 
ically based on magnitudes cuts (i.e., continuum domi¬ 
nated) and may miss faint galaxies with strong emission 
lines (i.e., large EWs). These biases can also affect the 
study of the SFR-M* relation since at low masses the suc¬ 
cess rate of getting redshifts may be higher if the SFR is 
above average whereas it may be lower for massive heav¬ 
ily attenuated galaxies (see, e.g., Rodighiero et al. 2014; 
Speagle et al. 2014, for a discussion of the biases and 
selection effects on the SFR-M* relation), thus introduc¬ 
ing a flattening in the SFR-M* relation. Medium- and 
narrow-band surveys are, by nature, less prone to this 
kind of selection biases. 

The so called ’’main sequence” (Noeske et al. 2007) 
of SF galaxies, has been extensively studied in the past 
using different sample selections from low- (see, e.g., 
Brinchmann et al. 2004; Elbaz et al. 2007; Ilbert et 
al. 2015) to high-redshift (see, e.g., Speagle et al. 2014; 
Steinhardt et al. 2014; Schreiber et al. 2015; Salmon et al. 
2015), in different environments, ranging from extremely 
low-density (known as cosmic voids; see Ricciardelli et 
al. 2014) to the most crowded regions of the Universe 
(i.e. filaments and galaxy clusters; see e.g., Koyama et 
al. 2013; Darvish et al. 2014; Vulcani et al. 2010, 2015), 
and using various diagnostics as SFR indicators (see, e.g., 
Rodighiero et al. 2014; Shivaei et al. 2015a,b). Despite 
the advancements on this topic, the debate on the nature 
and characteristics of the galaxies populating the main 
sequence is still open (see, e.g.. Kelson 2014, for a dif¬ 
ferent interpretation of the star-forming main sequence). 
As indicated above, the determination of the slope and 
the scatter of this relation and the comparison between 
different works is made difficult by various systematic 
effects (e.g., differences in calibrations, selection biases, 
used SFR indicator, among others; see, Speagle et al. 
2014 for a more complete discussion of these uncertain¬ 
ties). Furthermore, to disentangle between differences in 
the determination of the main sequence due to the men¬ 
tioned systematic effects, and pure evolutionary effects 
can be even more subtle. 

The main aims of this paper are: (1) to exploit 
the ultra-deep SHARDS data to select and investigate 
intermediate-redshift ELGs on the basis of their [OH] 
emission and (2) to perform an (almost) unbiased study 
of the SE activity of galaxies at intermediate redshifts 
(^1). This goal is attained by combining ultra-deep 
medium-band and multi-wavelength data (from UV to 
FIR) and comparing various SF indicators (UV, [OH] 
IR) to assess the robustness of the results. In particu¬ 
lar we focus on (3) the characterization of SFR-M* re¬ 
lation and (4) the dependence of the dust attenuation 
on the physical galaxy properties (M* and total SFR). 
The depth and quality of SHARDS data allow to im¬ 
prove upon previous intermediate redshift studies going 
down to fainter magnitudes and limiting the effect of the 
selection biases typical of spectroscopic surveys. 

The scheme of this paper follows. In Section 2, we 
summarise the main characteristics of the survey and the 
data used in this work. We introduce the sample selec¬ 
tion for [OH] galaxies and we define the different comple- 
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mentary samples (general stellar mass-selected sample, 
IR detected and UVJ/quiescent galaxies) used through¬ 
out the paper. We present the methods used to select 
ELGs, their redshift distribution and discuss the adopted 
procedures and caveats. In Section 3, we present the 
basic measurements for the equivalent widths and line 
fluxes characterizing the [Oil] population and discuss 
their observational properties. In Section 4, we discuss 
the [Ollj-based SFRs (SFR([OII])), the stellar mass-SFR 
relation, and their stellar population properties (mass 
and age distributions). In Section 5, we compare the 
properties for different sample selection and discuss the 
derived SFR-M* relations and the differences induced by 
each selection. Dust extinction and its dependence on 
the stellar mass and the SFR are investigated in Sec¬ 
tion 6. Main results and conclusions are summarised in 
Section 7. 

Throughout this paper we use AB magnitudes and the 
Ghabrier (2003) initial mass function (IMF). We adopt 
the cosmology Hq = 70 km s“^Mpc“^, 11^ = 0.3, and 
Da = 0.7. 

2. DATA AND SAMPLE SELECTION 

In this work, we have gathered a complete sample 
of star-forming galaxies at 2:~0.84 and z~1.23 in the 
GOODS-N field by combining 3 different selections: a 
pure stellar mass-selected sample of galaxies with SFRs 
calculated through the rest-frame UV emission at 150- 
250 nm, an emission-line selection based on the detec¬ 
tion of the [OII]A3727 line with data from the SHARDS 
(Perez-Gonzalez et al. 2013), and a selection of obscured 
SFGs based on mid- and far-IR data from Spitzer and 
Herschel. In the following subsections, we present the 
details about the different datasets used in the selection 
of SFGs and their characterization. We start with the 
description of the most unique selection among the three 
mentioned before: the [OH] sample selected with the new 
SHARDS data (see Sections 2.1-2.3). Then, we comple¬ 
ment this sample with the other selections and perform 
a comparison of the properties of the sample of [OH] 
emitters with the dust-obscured SFGs selected from IR 
data and the mass-selected sample (characterized by the 
UV emission, see Sections 2.4-2.5). This approach gives 
us a comprehensive view of the typical limitations inher¬ 
ent to the different selection methods and highlights the 
complementarity of the different selected samples. As 
we will discuss later, the construction of the complemen¬ 
tary samples must take into account the fact that the 
[OH] emitters selected with SHARDS have a very specific 
redshift distribution. Therefore, both the mass-selected 
and the obscured SFG samples should follow similar red¬ 
shift distributions. This is mandatory to be able to make 
a fair comparison between different samples and under¬ 
stand the SF and dust extincion properties of the whole 
population of SFGs at z~0.84 and Z'^1.23, which we will 
discuss in Sections 4, 5 and 6. 

2.1. Datasets gathered for this work 

GOODS-N is one of the most targeted areas of the sky 
at all wavelengths. For this work, apart from our new 
SHARDS data, we have combined the wealth of deep 
and high quality ancillary data, ranging from an ultra¬ 
deep X-ray exposure (2Ms GDFN, Alexander et al. 2003) 


to the deepest data in the MIR/FIR with surveys such 
as GOODS (Giavalisco et al. 2004), FIDEL (Frayer et 
al. 2006), PEP (Lutz et al. 2011), HerMES (Oliver et al. 
2010), and Herschel-GOODS (Elbaz et al. 2011). Multi¬ 
ple spectroscopic redshifts for faint targets are also avail¬ 
able (Wirth et al. 2004; Gowie et al. 2004; Reddy et al. 
2005, 2006; Barger et al. 2008; Gooper et al. 2011; Kriek 
et al. 2015). 

The SHARDS project was designed to be able to mea¬ 
sure absorption indices such as the Mg[/v , or D(4000) for 
galaxies at 2 = 1 .0-2.5 through imaging data, and detect 
ELGs up to 2~7. For those purposes, SHARDS obtained 
imaging data in the GOODS-N field through 25 medium- 
band filters covering the wavelength range between 500 
and 950 nm in sub-arcsec seeing conditions. The typi¬ 
cal width of our filters is 15 nm for wavelengths bluer 
than 880 nm, and 25-35 nm for the 3 reddest filters. In 
each filter, SHARDS is able to detect 26.5 mag sources at 
the 3 ct level (at least). Virtually all the deep region cov¬ 
ered by GOODS with HST/AGS is surveyed by SHARDS 
using two GTG/OSIRIS pointings, summing up a total 
surveyed area of ^130 arcmin^ (see Figure 2 from Perez- 
Gonzalez et al. 2013). The observations carried out by 
SHARDS allow to accurately determine the main prop¬ 
erties of the stellar populations present in these galax¬ 
ies through spectro-photometric data with a resolution 
R^50, sufficient to measure absorption indices such as 
the D4000 (e.g., Bruzual 1983; Balogh et al. 1999; Kauff- 
mann et al. 2003; Kriek et al. 2011; Hernan-Gaballero 
et al. 2013, 2014) or Mg[/v , index (Spinrad et al. 1997; 
McGarthy et al. 2004; Saracco et al. 2005; Daddi et al. 
2005; Cimatti et al. 2008; PCTez-Gonzalez et al. 2013). 
At this spectral resolution, it is also possible to detect 
emission lines and measure their fluxes and equivalent 
widths. Among the lines SHARDS can detect, we can 
mention Ha, [OH], [OHI], Lya, among others (see e.g. 
Villar et al. 2011; Perez-Gonzalez et al. 2013; Sobral et 
al. 2013). 

As explained in Perez-Gonzalez et al. (2013), given the 
special characteristics of the OSIRIS instrument at GTG, 
within each single frame taken with a given physical fil¬ 
ter, each pixel sees a different passband. To overcome 
this issue, we performed a detailed calibration as a func¬ 
tion of the position in the field of view (FOV). The sig¬ 
nificant variation of the passband seen by each point of 
the detector is a function of the position in the FOV and 
implies a complex behaviour of the absolute photometric 
calibration of the SHARDS images. To cope with these 
issues, we developed a special flux calibration procedure, 
aimed at determining the zeropoint of the SHARDS mo¬ 
saics in each filter as a function of position in the image. 
The flux calibration of the SHARDS mosaics was per¬ 
formed by comparing the measured photometry in our 
images with spectroscopic data for several sources in the 
field. 

To complement the SHARDS data, we also benefit- 
ted from the fact that GOODS-N has been observed 
by HST with AGS and WFG3 providing slitless, inter¬ 
mediate resolution spectroscopy in the optical (through 
the G800L grism; PEARS, Kiimmel et al. 2009; see also 
Pirzkal et al. 2004, 2009) and NIR (G102; PI: G. Barro, 
G141; PI: B. Weiner). In addition, the availability of the 
deepest IRAG ([3.6 /rm]<26.0 mag) and MIPS {F5a[2A 
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TABLE 1 

Characteristics of the SHARDS filters and observations used in 

THIS study 


Filter 

CWL 

Width 

Exposure time 

m 3 

a 

seeing 




PI 

P2 

PI 

P2 

PI 

P2 

( 1 ) 

( 2 ) 

(3) 

(4) 

(4) 

(5) 

(5) 

( 6 ) 

( 6 ) 

F670W17 

670.4 

16.0 

3795 

4554 

26.8 

26.9 

0.8 

1.0 

F687W17 

686.9 

17.2 

9270 

12360 

27.2 

27.1 

0.8 

0.9 

F704W17 

703.7 

17.9 

6120 

6120 

26.8 

26.8 

0.9 

0.9 

F806W17 

806.5 

16.1 

14900 

14900 

26.5 

26.6 

0.9 

1.0 

F823W17 

823.1 

14.7 

18570 

24760 

26.8 

26.8 

0.8 

0.9 

F840W17 

840.0 

15.6 

19530 

25872 

26.2 

26.4 

0.9 

0.9 


Note. — (1) Filter name. (2) Central wavelength (in nm) of the filter 
for angle of incidence AOI=10.5° (approximately that for the center of the 
FOV). (3) Filter width (in nm). (4) Exposure time (in seconds). (5) Average 
3-sigma depths (AB mag) for circular apertures of radius 0.8” for pointings 
1 and 2. ( 6 ) Average seeing (in arcsec) for pointings 1 and 2. 


^m]>30 /iJy) observations ensures the detection of the 
rest-frame NIR/MIR emission of the galaxies and allow 
us to estimate robust stellar masses and SFRs (Perez- 
Gonzalez et al. 2005, 2008; Cimatti et al. 2008). 

We merged all these dataset for a mass-selected sample 
as described in Perez-Gonzalez et al. (2008). In that pa¬ 
per, we built spectral energy distributions (SEDs) using 
aperture photometry and deconvolution algorithms for 
IRAC bands. We then fitted the SEDs with a set of tem¬ 
plates representative of the diverse galaxy populations 
and built with stellar population synthesis models. These 
fits allowed us to obtain estimations of the photometric 
redshifts and stellar masses. We describe in more detail 
this modeling in Sections 2.2 and 3. Based on the models 
best-fitting the photometric data, we also obtained lumi¬ 
nosities in the rest-frame UV (more specifically, at 150 
and 280 nm), which were then converted to SFRs using 
the recipes from Barro et al. (2011b), taken from Bell et 
al. (2005). Attenuation based on the UV slope /3 were 
also estimated using the Calzetti et al. (2000) attenua¬ 
tion law and the SFR(IR)/SFR(UV) vs. /3 calibration 
in Meurer et al. (1999). We will discuss in detail the 
UV-based SFRs in Section 5. 

2.2. Photometry and SED fitting 

In this section, we briefly summarise the basic steps 
and procedures adopted to extract the photometry and 
to fit the SEDs of different galaxy samples. These steps 
are performed exploiting the Rainbow Gosmological Sur¬ 
veys Database^ hosted by the Universidad Gomplutense 
of Madrid, which is the central repository for all the 
SHARDS and ancillary data used in this paper. 

Photometry is performed using the Rainbow (G. Barro 
et al. 2015, in preparation) tools that are able to carry 
out simultaneous photometry in multiple bands using 
prior based positions and apertures. The priors for the 
SHARDS extraction are based on SExtractor (Bertin & 

^ The Rainbow Cosmological Surveys database is a vast compila¬ 
tion of photometric and spectroscopic data for several of the deep¬ 
est cosmological fields, such as GOODS-North and South, COS¬ 
MOS, or the Extended Groth Strip, among others. It is publicly 
accessible through the website: https://rainbowx.fis.ucm.es/ 


Arnouts 1996) catalogs and Kron (1980) apertures ob¬ 
tained as an average of those measured in all SHARDS 
bands. Photometry for longer wavelength bands (start¬ 
ing with IRAG) is extracted using circular apertures of 
fixed sizes. 

For the SED fitting procedure, we adopt a two- 
population model, each population described by a SF 
history (SFH) following an exponentially declining law, 
characterized by timescales Tyou and a Tom parameter for 
the young and the old stellar population, respectively. 
We use the models from Bruzual & Gharlot (2003), with 
a Ghabrier (2003) IMF spanning stellar masses from 0.1 
to 100 Mq. We assume the dust attenuation law from 
Galzetti et al. (2000). For the whole sample of [OH] 
emitter candidates selected from SHARDS and the com¬ 
plementary samples described in Sections 2.3 and 2.4 we 
derive stellar masses, SFRs, dust reddening A(V), and 
stellar ages from the best fits. 

We compare the photometric data with a grid of mod¬ 
els probing r values in the range 6<log(T/2/r)<12 for 
each population. These r values were selected to include 
from instantaneous to roughly constant SFHs. The grid 
of models spans a range of ages, from 1 Myr to up to 
1 Gyr for the young population and from 1 Gyr up to 
the age of the Universe at the given galaxy redshift for 
the old population. We adopt a Galzetti et al. (2000) dust 
attenuation law and dust extinctions with Ayou(V)=0.0- 
7.0 mag and Aoid(V)=0.0-2.0 mag, with 0.1 mag incre¬ 
ments. We consider the six metallicities available from 
the Bruzual & Gharlot (2003) libraries. Additional de¬ 
tails and an example of the stellar population fits in the 
SHARDS spectral region are given in 3.1. In these fits, 
the redshifts are fixed to either the spectroscopic val¬ 
ues when available or the best photometric redshift solu¬ 
tion. As we show in Section 2.3.2, we have a considerable 
improvement on the photometric redshift determination 
using the full SHARDS dataset to perform SED fitting 
(G. Barro et al. 2015, in preparation; see also Ferreras 
et al. 2014). 

We remark that the results obtained from this model¬ 
ing are of particular relevance for the determination of 
the continuum at the wavelength of [OH] emission line, 
as discussed in Section 3. 
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Fig. 1. — Color-magnitude diagram showing ELG candidates at ~687 nm (left panel) and ~823 nm (right panel). The vertical axis shows 
the color between a central SHARDS filter band and the average magnitude in the adjacent (continuum) SHARDS bands. The locus for 
galaxies with an emission line within the filter width detected with more than 2.50- confidence is the region above the yellow-dashed curve, 
identified as the emitter locus. The vertical shaded orange line shows the minimum 3cr detection threshold of the SHARDS survey in the 
selected filter bands. Red dots indicate [OH] emitters candidates with spectroscopic confirmation, while blue diamonds stand for candidates 
selected using photometric redshifts. The red-dashed lines depict the typical photometric uncertainty as a function of the filter magnitude. 
Lower panels are a zoomed view around the selection curve. 


2.3. Sample selection of [Oil] emitters 

We have developed a novel selection technique to iden¬ 
tify ELGs using the SHARDS medium-band spectro- 
photometric dataset, based on similar narrow-band se¬ 
lection techniques. The technique to select star-forming 
galaxies (and also AGNs) from medium-band photome¬ 
try is based on comparing the flux measured in one filter 
(the central filter, hereafter) with those obtained in two 
adjacent filters (the continuum filters, hereafter). Galax¬ 
ies presenting emission lines would then pop-up in the 
central filter if their redshifts move the line within the 
passband of this filter, thus providing the emission line 
flux, while the adjacent filters would give us an estimate 
of the continuum around the emission line. 

A similar color-excess technique has been applied to 
multiple surveys of intermediate, high, and very high 
redshifts (focusing on different lines such as Ha or Lya) 
using typically narrow-band passbands as the central fil¬ 
ter and a broad-band passband as the continuum filter 
(see e.g. Ouchi et al. 2008; Villar et al. 2008, 2011; Sobral 
et al. 2009a,b, 2012, 2013; Matthee et al. 2014). 

Note that our dataset consists of medium-band filters, 


not narrow-band, with a typical equivalent spectral reso¬ 
lution i? ^ 50 to be compared, e.g., to i? 100 in Villar 
et al. (2008, 2011) or i? - 80 in Sobral et al. (2009b, 2012, 
2013), but the depth and image quality of the SHARDS 
images, jointly with the use of adjacent medium-band 
passbands for the determination of the continuum, offer 
several advantages over previous surveys as detailed next. 
Note also that the SHARDS spectral resolution is similar 
to that achieved with HST optical grism spectroscopy. 

2.3.1. Selection of ELG candidates using SHARDS data 

In this work, we used two SHARDS filters as central 
passbands, F687W17 and F823W17, two of the deep¬ 
est in our survey, centered at 687 and 823 nm, respec¬ 
tively. The most common emission line detected with 
these filters is [OHJA3727, which means that the central 
wavelength bf corresponds approximatively to a redshift 
of z^O.84 and 2^1.23 for the bulk of our emission-line 
galaxy sample. We start the analysis of SHARDS se¬ 
lected ELGs with these two filters because (1) they are 
representative enough of the whole dataset, (2) they are 
two of the deepest ones, and (3) they enable the study of 
samples of galaxies located at sufficiently separated red- 
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Fig. 2.— Spectroscopic and photometric redshift distributions for the sources selected as ELGs candidates with the F687W17 and 
F823W17 filters. The histograms in red refer to galaxies with spectroscopic redshift confirmation. The outlined black histograms refer to 
photometric redshifts (Perez-Gonzalez et al. 2008; G. Barro et al. in preparation). The expected redshifts for emitters with some of the 
most typical lines (e.g., Lya, [Oil], [OIII], or Ha) lying in each one of the two central filters are marked with shaded regions according to 
the legend of each panel. Numbers below the histogram distributions help to recognise each line. We also mark other spectral features such 
as the Mg-Fe absorption band at ~280 nm, or the 2175 A dust absorption bump, which may mimic emission lines at a different redshift. 
Bin sizes are proportional to A( 2 )/(l + 2)~0.006, our typical photometric redshift accuracy. 


shifts in order to unveil the possible evolutionary effects 
on their derived physical properties. In a future work 
we will present more extensive results including all the 
SHARDS filters. 

In order to be able to compute the continuum level, 
we have also used data in four additional filters (lying on 
the bluer and redder side of each central filter), namely: 
F670W17, F704WI7, F806W17, and F840WI7. The 
main characteristics of the central and continuum filters 
and observed data used in this work are summarized in 
Table I. 

A sketch of the technique used in the selection of the 
ELGs sample is presented in Figure I. We compare the 
magnitude in the central filter (F687WI7 or F823WI7) 
with the continuum filters (a combination of the 2 ad¬ 
jacent ones for each central filter). Typically, the esti¬ 
mation of the continuum in ELGs is carried out with 
broad-band filters, which also includes the line entering 
in the central filter and possibly others Villar et al. (2008, 
2011); Sobral et al. (2009b, 2012, 2013). By using two 
adjacent medium-band filters, we can estimate a more 
accurate continuum level. This method can be consid¬ 
ered robust for the continuum determination as far as the 
contribution of the emission line to the side filters can be 
neglected, which is typically the case for SHARDS filters 
(given the filter shape). Indeed, the SHARDS contin¬ 
uum coverage of the optical spectral range ensures that 
intrinsic variations of the spectral slope do not affect the 
continuum determination. The probability of another 
line entering in the continuum filter is also smaller than 
when using broad-band filters. Nonetheless, the finite 
width of the filters, the noise, and other spectral features 
possibly falling in our observed spectral range can affect 
the final measurement. We rely on the interpolation be¬ 
tween the two filters to define the observed continuum 
for the selection of ELGs. However, once the emission 
line candidates are selected, we perform a detailed SED 
fitting to fine-tune the continuum determination in or¬ 
der to obtain robust line fluxes and equivalent widths 
(as described in Sections 2.2 and 3.1). 


Figure 1 shows the selection diagrams for the filters 
F687W17 (left panel) and F823W17 (right panel). They 
depict the ’’continuum—line color” vs. the ”central-filter 
magnitude” for each of the two filters. As we move 
to fainter magnitudes, photometric errors become larger 
and larger, resulting on a typical “trumpet” shape of the 
cloud of points. The color vs. magnitude distribution 
can be used to estimate the typical scatter as a function 
of magnitude. In this plot, a is the standard deviation 
of the colors, and we calculate it as a function of magni¬ 
tude. Galaxies with significantly larger colors than the 
statistical standard deviation (associated to the photo¬ 
metric uncertainties) are identified as sources with an ex¬ 
cess of light in the central filter, i.e., emission-line galaxy 
candidates. We have chosen a selection function (yellow 
curve) coming from a 2.5 ct cut in the color-magnitude dis¬ 
tribution. The vertical orange line represents the mini¬ 
mum detection threshold adopted for the selection, corre¬ 
sponding to the minimum 3a detection limit of the three 
filters used in each case (selection band plus the two con¬ 
tiguous filters, see Table 1). These cuts, provide a good 
compromise between the inclusion of the largest number 
of ELG candidates and the contamination from spurious 
sources (see, e.g., Villar et al. 2008, 2011). Summariz¬ 
ing, galaxies identified as potential emitters are found in 
the emitters locus region of these diagrams. Note that 
the selected galaxies are candidates to have an emission 
line located in the central filter. This line could be, for 
the F687W17 filter, [OH] at z^O.84, but it may also be 
any other line such as Lya at z^6 or Ha at z^O, among 
others. 

The redshift distributions of the galaxies found in the 
emitter loci presented in Figure 1 are shown in Figure 2. 
Most of the peaks in the distribution are actually due to 
an emission line falling within the central filter observa¬ 
tional window. Spikes can also be produced by absorp¬ 
tion lines and/or breaks in the continuum falling in one 
of the filters used for the continuum definition, which can 
mimic an emission line at a different redshift (see Fig. 9 
in P&ez-Gonzalez et al. 2013). For example, a number 
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Fig. 3.— Comparison between spectroscopic and photometric 
redshifts for SHARDS sources around 2:'^0.84 {top panel) and 
Z'^1.23 {bottom panel), the two redshifts where we present [Oil] 
emitters selected with the F687W17 and F823W17 filters. We show 
the distributions of redshift accuracies, = {Zpj^o^ — Zspec)/{^ + 
Zspec)- In both panels, the blue curve is the Gaussian fit to the 
distribution, with the vertical solid and dashed lines representing 
the mean offset and Icr dispersion respectively. 


of galaxies can be selected due to the D4000 break (see 
Hernan-Caballero et al. 2013) falling in the filter used 
to define the left side of the continuum, mimicking an 
emission line at lower redshift. 

The photometric redshift distribution broadly follows 
the shape outlined by the spectroscopic sample, reveal¬ 
ing their high quality. In Figure 3, we show a comparison 
between the SHARDS-derived photometric redshifts and 
available spectroscopic redshifts in the GOODS-N field 
for our two redshift ranges around 2~0.84 and 2~1.23. 
We have estimated a typical accuracy of ^{z)/{l -|- 
z)=0.006 by comparison of the photometric redshift with 


their corresponding spectroscopic measurements at all 
redshifts and magnitudes. In the intervals covered by 
the two filters used in this work, we derived a typical 
rms ~0.5%. The high quality of the photometric red¬ 
shifts is very relevant when aiming at selecting a partic¬ 
ular class of ELGs, such as [Oil] emitters at two different 
redshifts, in our case. These photometric redshifts will 
be presented in another work (G. Barro et al. 2015, in 
preparation). 

For this study, we focus only on the [Oil] emitters se¬ 
lected combining Figures 1 and 2. To make the most 
robust selection, the tolerance allowed for the inclusion 
of [Oil] ELG candidates is different depending on the 
type of redshift estimate at hand. Galaxies with already 
available spectroscopy are easily confirmed as [Oil] emit¬ 
ters selected by SHARDS. The selection of galaxies with 
no spectroscopic confirmation must be done with pho¬ 
tometric redshifts, accounting for their uncertainty (see 
the next section for more details). 

2.3.2. Identification of [Oil] emitters candidates 

Therefore, the question is how we can identify [OH] 
emitters and separate them from galaxies with other 
emission lines or bumps. We identify bona fide [OH] 
emitters as those ELG candidates for which the spec¬ 
troscopic redshift is such that the [OH] line lies within 
the interval implicitly defined by the FWHM of the ob¬ 
served central passband. Note that, as explained in Sec¬ 
tion 2.3 and Perez-Gonzalez et al. (2013), the actual pass- 
band seen by each galaxy depends on the exact position 
in the GTG/OSIRIS FOV. The procedure is graphically 
explained in Figure 4. The medium-band technique is 
very effective in selecting ELGs in a narrow range of 
redshifts. If we had spectroscopic confirmation for all 
sources, those ranges would be as narrow as A(2)~0.05. 
However, ^ 53% of of galaxies finally identified as [OH] 
emitters at z^O.84 and ^ 74% of the sample z^l.23, 
only count on a photometric redshift. In order to de¬ 
termine which galaxies are actual [OH] emitters on the 
basis of photometric redshifts, we have to consider the 
error of those photo-z’s. 

In addition to this, due to the variable central wave¬ 
length (CWL) along the FOV of the OSIRIS instrument, 
we obtain a continuum distribution for the selection func¬ 
tion that depends on the CWL of the filter for each 
galaxy position. Combining all effects, we obtain the 
hatched region in Figure 4, which represents the overall 
extension in the redshift space for our selection function 
based on the trumpet plots in Figure 1 and the photomet¬ 
ric redshifts (see Figure 2). We remark that this redshift 
range does not apply to the selection in the whole area 
surveyed by SHARDS, the redshift range probed actually 
depends on the position of the galaxy in the FOV. 

This procedure provides two catalogues of ELGs se¬ 
lected as [OH] emitters based on the two redshift es¬ 
timates: one using photometric redshifts, and one on 
the basis of spectroscopic redshifts, which should be a 
subsample of the former. We call the latter the specz- 
confirmed ELGs, while the former is the photoz-selected 
sample. These catalogues can be merged to get a robust 
and complete selection. Adopting a similar procedure, 
we are able to produce catalogs for each possible detected 
line (as those listed in the legend of Fig. 2) and we plan 
to analyze the full combined set of detected lines in a 
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Fig. 4.— Scheme of the redshift dependence of the selection of [Oil] ELGs based on the SHARDS medium-band data. We show the 
results for filters F687W17 (2:'^0.84) and F823W17 (2:'^1.23). Blue (red) dashed lines show the transmission curves for the bluest (reddest) 
passbands seen by any galaxy in the SHARDS images for these filters. The orange curves represent the nominal central passband for the 
given SHARDS filter (centered on 687 and 823 nm). The filled grey area represents the interval where galaxies with spectroscopic redshifts 
are selected. The shaded pink area indicates the redshift interval of the selection for galaxies with photometric redshifts. Its width is based 
on the photometric redshift mean error < 5(z) >= cr x (1 -|- 2 ) at redshift 2 , where a = 0.005 for 2^0.84 and a = 0.006 for 2:'^1.23. The 
50% normalized transmission is marked with a black dashed lines. The grey-hatched area gives the global envelope of the transmission 
curves along the FOV. 


future work. 

In order to assess and improve the selection of [Oil] 
emitters, we compare the spectroscopic and photomet¬ 
ric selected samples in Figure 5. In the top panels, we 
plot the redshift distribution of the two samples of [Oil] 
emitters based on spectroscopic and photometric red¬ 
shifts. The photometric redshift selection closely follows 
the spectroscopic distribution, but it includes a higher 
number of galaxies. If the galaxies selected with photo¬ 
metric redshifts are bona fide [Oil] emitters, then half of 
the 2^0.84 sample and one fourth of the 2^1.23 sample 
currently would have a spectroscopic confirmation (i.e., 
the spectroscopic completeness is ^50% and ~25%). 

It is interesting to note that the galaxies with no spec¬ 
troscopic confirmation are typically fainter than AB~ 
24.5, as clearly shown in Figure 1. In fact, the spec¬ 
troscopic completeness for the selected bona fide [Oil] 
emitters brighter than R = 24.5 is 91% (74%) in the 
F687W17 (F823W17) filter, and ~ 13% (16%) for galax¬ 
ies with i?>24.5 mag. This is the spectroscopic limit of 
the redshift surveys carried out in the GOODS-N held, 
and the typical detection threshold for the vast major¬ 
ity of data taken with state-of-the art spectrographs in 
10-meter class (or smaller) telescopes. The SHARDS ob¬ 
servations reach at least 2 magnitudes fainter, and thus 
open the possibility to reliably select and study fainter 
and/or higher redshift ELGs. We note, furthermore, that 
the number of galaxies having publicly available spectra 
of good quality, where the S/N is sufficiently high to de¬ 
rive a robust equivalent width (EW) and [Oil] line flux 
determination, only represents a small fraction (typically 
^ 5 — 10%) of the total number of galaxies for which a 
spectroscopic redshift determination is available in this 
study. 


In the second and third rows of Figure 5, we compare 
more quantitatively the samples selected with spectro¬ 
scopic and photometric redshifts in terms of the success 
rate and the contamination. The success rate is defined 
as: 

SR = Nph-conf/Nsp (1) 

where Ngp is the number of galaxies in the spectroscopic 
catalog and Nph-conf is number of galaxies with a spec¬ 
troscopic confirmation which are included in the photo¬ 
metric redshift sample. This number is related to the 
ability to select a galaxy using the photometric redshift 
and provides an estimate of the number of galaxies we 
lose due to uncertainties in the photo-z or a failure in 
the detection of the emission-line with the SHARDS fil¬ 
ter. Typically, the selection based on SHARDS medium- 
band data and photometric redshifts has a success rate 
higher than 90%. 

In an analogous way, we define the contamination as: 

G = 1 Mpfi—cQjif/JVph—sp (2) 

where Nph-sp is the total number of galaxies in the pho¬ 
tometric catalog having a spectroscopic redshift estimate 
(either within the expected redshift range or outside). 
This number gives an estimate of the fraction of contam¬ 
inants expected in the final sample. The contamination 
measured in both redshifts bins considered in our study 
is less than 10%. The measured success rate and con¬ 
tamination level are coherent with the expectations due 
to the photometric redshift errors and spectroscopic un¬ 
certainties for the faintest observed galaxies. 

In the bottom panels of Figure 5, we show the red¬ 
shift distribution for the final sample of bona fide [OH] 
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Fig. 5.— Redshift histograms showing the selection, success rate and contamination rate for the [Oil] ELG candidate samples built with 
central filters F687W17 (2:~0.84, plots on the left) and F823W17 (-3^1.23, right plots) using spectroscopic or photometric redshifts. Top 
panels: redshift distributions for the selected objects using spectroscopic (red filled histogram) or photometric (black histogram) redshifts. 
Second row of panels: redshift distribution of objects selected on the basis of their spectroscopic redshift (red filled histogram) compared 
to the photometric redshift distribution for the same galaxies (black histogram). The ratio in the numbers of objects in these distributions 
defines the success rate (SR). Third row of panels: redshift distribution for objects selected on the basis of their photometric redshift (black 
histogram) compared to the distribution of their corresponding spectroscopic redshifts (red filled histogram) for those objects in which 
both redshift estimates are available. The normalised difference in the numbers of objects in these distributions defines the contamination 
(C). Bottom panels: redshift distribution for the final selected sample built by complementing the spectroscopic and photometric redshift 
selections. The vertical dashed line represents the median redshift in each panel. 
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Fig. 6 .— Top panels: UVJ color-color diagram with the quiescent galaxies region defined as in Whitaker et al. 2011 (F687W17 to the left, 
F823W17 to the right). IR detected galaxies within the UVJ quiescent region are excluded from the passive sample. The different samples 
are indicated with different symbols identified in the legend. Bottom panels: sSFR-M« diagrams for the two SHARDS filter selections at 
z~ 0.84 (left panel) and z~ 1.23 (right panel). The different subsamples are plotted with different symbols identified in the legend. UVJ 
passive galaxies are defined as in Whitaker et al. 2011. The dark-grey dashed and dotted horizontal lines mark two cuts at constant sSFR 
(log(sSFR [yr“^])=—9.7 and —9.5), adopted for the sSFR-based definition of quiescent galaxies. The black-dotted lines mark the sSFR-M* 
relations for constant values of SFRs (SFR=3 and 6 Maun yr~'^)i corresponding approximatively to our detection limits for IR-selected 
galaxies in the two redshift ranges. 


emitters built by merging the samples selected on the 
basis of spectroscopic and photometric redshifts. These 
final catalogs do not include the known contaminant 
galaxies identified on the basis of measured spectro¬ 
scopic redshifts. Accounting for this, the final contami¬ 
nation fractions are 2% and 6% for the filters F687W17 
and F823W17, respectively. Obviously, both the con¬ 
tamination and the success rate for the final selection 
are very sensitive to the photometric redshift errors. 
We have verified the effect of using photometric red¬ 
shift catalogs derived using only broad-band photome¬ 
try (that is, excluding SHARDS data), with an accuracy 
of A(z)/(1 + z)'^0.03, typical of most photometric red¬ 
shift catalogs. The success rate decreases to 70% level, 
and the contamination increases to ^ 40%. Henceforth, 
we will always refer to the sample shown in the bottom 
panels of Figure 5 as the final [OH] emitters sample. 

2.4. Complementary samples: mass-selected, 
IR-detected and quiescent galaxies 


As already mentioned before, we use in our analysis the 
ancillary multi-wavelength catalog and advanced prod¬ 
ucts in the GOODS-N field presented in Perez-Gonzalez 
et al. (2008) and compiled in the Rainbow Cosmologi¬ 
cal Surveys Database (see, Perez-Gonzalez et al. 2008; 
Barro et al. 2011a). This dataset includes observations 
from X-rays to the Far-IR and radio bands, as well as 
spectroscopic data in the GOODS-N field from the lit¬ 
erature. In Perez-Gonzalez et al. (2008), we presented a 
merged photometric catalog including broad-band data 
for a stellar mass selected sample based on ultra-deep 
IRAC observations. For this work, we have merged this 
catalog with the SHARD dataset. As described in Perez- 
Gonzalez et al. (2008), the parent sample used in this pa¬ 
per has been built using IRAC luminosities as proxies for 
the stellar mass. Thanks to the ultra-deep IRAC imaging 
available our sample is mass-complete for galaxies with 
M>109 Ms„„ up to redshift z=l. 

In order to perform a comprehensive study of SFGs at 
z~0.84 and z~1.23 , we have built a general parent sam¬ 
ple, starting from the GOODS-N mass-selected sample, 
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by selecting all the galaxies that would enter the medium- 
band hlter selection using the SHARDS data (i.e. with 
the same redshift distribution as our [OH] emitters), in¬ 
dependently from their emission properties. 

We apply the same selection procedure presented in 
Section 2.3.2 for the selection of [011] emitters, but this 
time relaxing the condition on the emission properties 
of the galaxies. This means that the position-dependent 
redshift selection function (as schematically represented 
in Figure 4) will not be applied only to the ELG candi¬ 
dates (i.e. galaxies within the emitter locus in Figure 1) 
but to the whole mass-selected galaxy sample, indepen¬ 
dently from the emission properties. In other words, we 
check the position within the FOV of each galaxy se¬ 
lected from its (photometric or spectroscopic) redshift, 
and then determine whether that redshift is good for a 
selection based on the SHARDS filters should the galaxy 
had had an emission line. In this way we can define a 
parent sample that includes and complements our [OH] 
emitter compilation and provides a complete census of 
galaxies (both star-forming and also quiescent) at those 
redshifts, which we simply dub the parent sample here¬ 
after. 

Starting from this general selection, we can derive dif¬ 
ferent subsamples according to several indicators for the 
SFR. In fact, we can estimate the SFR from the rest- 
frame UV luminosity (corrected for attenuation) derived 
from the SED fitting for all the parent sample, then we 
can define the subsamples of [OH]-detectecl (according 
to our SHARDS selection) and IR-detected galaxies (for 
example, using Spitzer/MIPS 24 ^m data). Finally, we 
can define a sample of quiescent galaxies using the rest- 
frame (V—J) versus (U—V) color-color diagrams follow¬ 
ing Whitaker et al. (2011), or alternatively as all galaxies 
lying below a certain sSFR cut (we use 0.2 Gyr“^ as our 
arbitrary limit) and not presenting IR emission. 

Interestingly, the sample of quiescent galaxies defined 
on the basis of their sSFR does not fully overlaps with 
the definition of UVJ-passive galaxies usually adopted 
in literature (Whitaker et al. 2011), although relaxing 
the cut in sSFR would improve the agreement, espe¬ 
cially for the higher redshift sample. In fact, increas¬ 
ing the cut from log(sSFR [yr“^])=—9.7 to log(sSFR 
[yr“^])=—9.5 (shown as horizontal dashed and dotted 
dark-grey lines in the bottom panels of Figure 6 respec¬ 
tively), the fraction of commonly selected objects would 
rise from — 55% to — 85% for the z — 1.23 sample. Here¬ 
after, we will refer as the ’’quiescent sample” to the one 
selected on the basis of the UVJ diagram, unless stated 
otherwise. Summarizing, the number of galaxies in each 
sample is 585/285/143/47 and 332/142/83/21 for the 
comparison/[OH]/IR/UVJ subsamples in the F687W17 
and F823W17 filters respectively. The definition of the 
samples is illustrated in Figure 6. We observe also 
that our IR detection limit is SFR—3 Mgun yr~^ and 
SFR-6 M,u„ yr-i, for the F687W17 and F823W17 fil¬ 
ters respectively. We mark these constant SFR relation 
as black-dottet lines in Figure 6. 

In Fig. 7, we show the spectroscopic and photometric 
redshifts distributions for the total population of galaxies 
in GOODS-N extracted from the Rainbow Database. We 
compare this global distribution with the redshift distri¬ 
bution of bona fide [OH] emitters identified through the 



REDSHIFT 

Fig. 7. — Histograms illustrating the redshift distributions for 
the selected sample of the [Oil] emitters (blue shaded histogram) 
and the parent sample (yellow shaded histogram) extracted from 
the parent catalog applying the same filter selection function. For 
completeness also the photometric and spectroscopic redshift dis¬ 
tributions for the parent catalog are shown (black/red outlined 
histograms respectively). Note that for the lower redshift selection 
(upper panel) the presence of a prominent large scale structure is 
evident. 

two SHARDS filters used in this work (see Section 2.3.2) 
and with the parent sample defined above. From this 
figure we notice that in the case of filter F687W17, cor¬ 
responding to the redshift slice at 2 — 0.84, we are se¬ 
lecting our [OH] emitters from a large scale overdensed 
region, clearly identified by the global spectroscopic (and 
also photometric) redshift distribution. Note also that 
the photometric and spectroscopic redshift distributions 
closely agree, confirming the goodness of photometric 
redshift estimates. 

2.5. AGNs fraction 

Our method to detect ELGs is sensitive to both SFGs 
and AGNs. Since we are interested in the former, we 
discuss here the possible contamination from AGNs in 
our sample. Of course, these AGNs likely have also active 
SF. 

The SHARDS field fully overlaps with the 2 Ms Ghan- 
dra Deep Field North (CDFN; Alexander et al. 2003). 
The X-ray source catalog from Alexander et al. (2003) 
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contains 267 X-ray sources in the area surveyed by 
SHARDS. The Chandra images are deep enough to de¬ 
tect the X-rays associated with normal SF (e.g., Per¬ 
sic et al. 2004; Lehmer et al. 2010; Pereira-Santaella et 
al. 2011). In order to avoid contamination from non- 
AGNs X-ray sources, we restrict our AGN identification 
to sources with X-ray luminosities Lx > 10^^ erg s“^ 
from Barger et al. (2008) in either the soft (0.5-2 keV) 
or hard (2-8 keV) band. Cross-matching our samples of 
[Oil] emitters with the ultra-deep X-ray CDFN catalog, 
we identify 4 AGNs in the F687W17 sample at 2^0.84 
and 2 AGNs in the F823W17 sample at 2:~1.23. These 
numbers correspond to a fraction of '^1.2% of the sam¬ 
ples, values in good agreement with the X-ray selected 
AGN fraction of 1 — 2% estimated by Zhu et al. (2009) in 
their analysis of z^l [Oil] emitters in the DEEP2 sur¬ 
vey. Giardullo et al. (2013) found a slightly higher frac¬ 
tion of AGNs (~ 3%) from their sample of [011]-selected 
galaxies in the HETDEX pilot survey (Hill et al. 2008; 
Adams et al. 2011), although they used a fainter limit 
for identifying AGNs from the parent X-ray catalog. As 
mentioned above, the correlation between X-ray emis¬ 
sion and AGN contribution to the derived SFR does not 
necessarily imply that the integrated [OH] emission is 
dominated by flux from the central engine. However, the 
SFR estimations based on indicators such as the [OH] 
emission could be contaminated by the AGN. Indeed, 
by examining the few X-ray bright objects in our survey 
we verified that these candidate AGNs are generally not 
amongst the strongest [011] emitters. The uncertainty 
of contribution of the AGN to the observed [OH] emis¬ 
sion and the small number of objects involved lead us 
to conclude that the exclusion of candidate AGNs from 
the final sample constitute a negligible source of uncer¬ 
tainty in our results. We note that a small fraction of 
highly obscured AGN could be missed by selecting only 
on the basis of the observed X-ray emission, which is 
able to select mostly unobscured and moderately ob¬ 
scured AGN (see, e.g., Fiore et al. 2009, 2012). How¬ 
ever, given the extreme depth of the available X-ray data 
in the GOODS-N field, and the already modest fraction 
(^ 1%) of unobscured-l-moderately obscured AGN de¬ 
tected in our sample, we can assume that the additional 
number of highly obscured AGN is negligible. 

3. OBSERVATIONAL PROPERTIES OF THE [OH] 
EMITTER SAMPLE 

In this section, we describe the procedure to measure 
[OH] fluxes and EWs from the SHARDS data for the 
[OH] sample of galaxies described in the previous section. 
The method is based on the estimation of the contin¬ 
uum emission around the emission line by comparing the 
SHARDS spectro-photometric data with stellar popula¬ 
tion models, and then using the central filter to measure 
the line flux. This modeling is described in Section 2.2. 
Then we compare the properties of our [OH] emitters 
with other similar samples from the literature. 

3.1. Line fluxes and equivalent widths o/[OH] emitters 

For the final samples of [OH] emitters at 2^0.84 and 
z^l.23 selected from the F687W17 and F823W17 filters, 
we have measured line fluxes and EW using the SHARDS 
medium-band data as detailed next. 


The flux emitted in the [OH] line is by definition (see, 
e.g., Villar et al. 2008; Sobral et al. 2012) : 


F([OH]) = (/^ 


central 


jcontinuum 


) X A 


(3) 


where /“ntraj jcontinuum measured flux den¬ 

sities in the central filter containing the emission line and 
the continuum obtained from the stellar population fits 
(defined as described below in this section). is 

the FWHM of the SHARDS central filter. The derived 
observed equivalent width is: 


EW{[Oll]) = ^([011])//^°"*“““™ (4) 

where F([OH]) is the line flux defined above. 

We have followed a careful and sophisticated method 
to measure EWs and fluxes, based on a robust determi¬ 
nation of the continuum around the emission line. The 
continuum estimation can be extracted from a linear in¬ 
terpolation between the two continuum hlters used in the 
selection of ELG candidates, but this method can be af¬ 
fected by line contamination of the side filters or other 
local features. A more accurate estimation relies on the 
use of stellar population models htting the medium-band 
data around the emission line (see, e.g., ?Vilella-Rojo et 
al. 2015, for a similar approach), thus exploiting all the 
information available from the SHARDS data. 

In Figure 8, we illustrate the adopted procedure for 
a galaxy at z~0.84. The SHARDS data are fitted to 
stellar population synthesis models as described in Sec¬ 
tion 2.2. When performing the fit, we excluded the filter 
containing the emission line, which is obviously contami¬ 
nated. The model allows us to determine the continuum 
at the position of the emission line. The continuum ob¬ 
tained from the best fit (green segment in Figure 8) is 
usually lower (on average ~10-15% lower density flux) 
than that obtained from the linear interpolation (grey 
line) of the two adjacent hlters. This translates into sys¬ 
tematically larger (but more robust) EW for the line (on 
average '^35-40% larger EWs) and higher huxes (on av¬ 
erage ~25-30% larger huxes). In the lower panel of Fig¬ 
ure 8, we show the effect of including or excluding the 
hlter that contains the emission line. Typically, there 
is no much diherence between the two hts (< 10% dif¬ 
ference in derived EW), but larger discrepancies (up to 
50% difference in derived EW) can be found when 
the emission line is very strong, affecting the continuum 
determination. 

In order to assess the reliability of our measurements 
based on spectro-photometric data from SHARDS, we 
compare line huxes derived from our medium-band hl¬ 
ters to those derived from spectra of the same objects 
and available from literature (TKRS, Wirth et al. 2004, 
and DEEP3 data. Cooper et al. 2011). The subsample 
of galaxies having spectra of good quality, i.e. S/N high 
enough to robustly measure EW([OH]) and F([OH]), 
contains 42 galaxies (~10% of the hnal sample). We 
show the comparison between the line huxes measured 
on the basis of the publicly available spectroscopic and 
SHARDS photometric data in Figure 9. 

On average, huxes measured from spectroscopy are 
~30-35% larger than huxes measured from SHARDS 
photometry when using a continuum determined from 
linear interpolation between adjacent hlters. Our proce¬ 
dure to estimate the continuum based on stellar popula- 
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Fig. 8. — Illustration of the method used to refine flux and EW 
measurements from SHARDS data by determining an accurate con¬ 
tinuum level (upper panel) for [Oil]. For all the galaxies selected 
as [Oil] emitters, we fit the SHARDS medium-band data to stel¬ 
lar population models and extract the continuum level from the 
best fitting model. Black filled stars are the SHARDS medium- 
band observed fluxes (with the associated 1—cr error bars), with 
the larger blue star indicating the filter used for the detection of 
the [OH] emission line. In red, we show the best-fit model to the 
SHARDS data (see text for details). Open stars represent the 
flux predicted from the best-fit model, after convolution with the 
SHARDS filters. We mark with a horizontal grey segment the 
continuum obtained by simple interpolation of the two adjacent 
filters to the central passband. The continuum obtained from the 
model best-fit is marked in green. We also illustrate (lower panel) 
the effect of excluding (continuous red line, same as in the upper 
panel) or including (dashed blue line) the photometric data-point 
corresponding to the filter containing the [OH] line. The small 
inset shows the global two-population fit (in red) depicted in the 
two main panels, with the young/old components overplotted in 
blue/yellow respectively. 


tion synthesis models significantly improves the compar¬ 
ison between spectroscopic and photometric line fluxes, 


Fig. 9. — Comparison between spectroscopic- and SHARDS- 
derived line fluxes and EWs for the [OH] line (upper and lower 
panel respectively). Spectroscopic measurements are obtained from 
public data (TKRS and DEEPS spectra). The continuum grey line 
represents the identity relation. Blue points and orange diamonds 
represent flux and EW estimates based on best-fit or linear in¬ 
terpolation for the continuum, respectively. The blue-dashed and 
orange-dotted lines are the linear best-fit to the two datasets. 


making them consistent within ^5-10%. Despite the 
large scatter, this result supports the robustness of our 
approach. Note that the spectroscopic measurement of 
the line flux is affected by slit losses, and typically this 
correction is carried out by assuming that the whole 
galaxy has a constant EW equal to the one measured 
in the spectra. If galaxies count with nuclear bursts, this 
correction would overestimate the total flux of the galaxy. 
This is confirmed by Figure 9, where we show that the 
observed EW([OII]) from SHARDS data are, on aver¬ 
age, ^30% smaller than those measured in spectra. In 
contrast, photometric measnrements of the line flux us- 
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ing integrated magnitudes are not affected by aperture 
effects, although spatial gradients of the emission could 
affect the selection of ELG candidates. 

We show the resulting distribution of rest-frame 
EWs (EWrf([OII])), F([OII]) and i?-band magnitudes 
for the final sample in Figure 10. Errors in the es¬ 
timates of the line flux and continuum level propa¬ 
gate into uncertainties in the derived EWrf([OII]) of 
^^20% at most. We observe a median (and [quar- 
tiles]) rest-frame EWrf ([01I])=51[38,65] A at 2^0.84 and 
EWrf([OII])=67[49,91] A at z^l.23. A -0.5 mag shift is 
also seen in the observed i?-band magnitude distribution, 
with the higher redshift sample peaking at fainter mag¬ 
nitudes. The line flux distribution is similar in the two 
cases, ensuring that we are not introducing any partic¬ 
ular redshift-dependent bias in line flux measurements. 
The median observed line flux for both distributions is 
—2x10“^^ erg/s/cm^, with measurements extending to 
—2x10“^® erg/s/cm^ on the faint end. We note that the 
minimum rest-frame equivalent width achievable assum¬ 
ing the typical uncertainty for the faintest [011]-selected 
galaxies corresponds to 30 A and 36 A for the higher and 
lower redshift sample respectively. For brighter galax¬ 
ies with smaller uncertainties we can reach even smaller 
values for the EW. The smaller value for the higher 
redshift sample is due to the smaller rest-frame filter 
width at this redshift. These values are indicated as a 
dot-dashed vertical blue line in the left panels of Fig¬ 
ure 10. These lower limits ensure a robust determination 
of the median EWrf ([Oil]) values discussed in the fol¬ 
lowing section. Indeed, we can convert these values to 
a lower limit in the line-flux detection (F([Oil])limit^^6- 
8xl0“^®ergs“^ cm“^), also reported as vertical blue 
dot-dashed line in the central panels of Figure 10. In 
this case, as expected, a slightly higher observed line flux 
limit is found for the higher redshift sample. 

These flux and EW detection limits are comparable 
to those obtained with deep spectroscopy. We conclude 
that our SHARDS spectro-photometric survey is very ef¬ 
fective in selecting ELGs down to the faintest flux levels 
achieved by spectroscopy. Note that the median fluxes 
and EWs of our samples are significantly higher than the 
detection limits, so the results presented in the following 
sections are not significantly affected by incompleteness 
at the faintest flux or smallest EW levels. 

3.2. EWrf([OII]) evolution 

The measurement of the EW([OII]) provides a valu¬ 
able source of information on the relevance of the ongo¬ 
ing star formation in galaxies. As it is well known (e.g., 
Kennicutt 1998) the [Oil] line flux is strictly related to 
the radiation field coming from the young {t < 20 Myr) 
and massive stellar (M> 8 Mgun) populations (and also 
to metallicity and other properties such as density and 
temperature). On the other hand, the continuum is more 
representative of a larger timescale (t—1 Gyr) star for¬ 
mation. The ratio of these quantities, i.e. the EW([OII]), 
thus provides an estimate of the current SF efficiency at 
the epoch of observations with respect to the average SF 
in the past galaxy life. This is an indicator of the current 
galaxy SF efficiency and it is largely independent on the 
degree of internal extinction of a galaxy. 

In some recent works (see, e.g., Ciardullo et al. 2013; 


Lee et al. 2012; Kornei et al. 2012; Ly et al. 2012) the 
distribution of the EW for large samples of [Oil] ELGs 
are presented and their median values given for samples 
ranging from z=0 to z—I. Locally, the distribution of 
rest-frame EW([OII]) for [OII]-selected galaxies is known 
to peak near 5 A, and then it rapidly decays without any 
noticeable dependence on the galaxy luminosity (Blanton 
& Lin 2000). 

In Figure 11 we show the evolution of the me¬ 
dian EWrf([OII]) measurements of [OII]-selected sam¬ 
ples from literature at redshifts 0<2:<1.5 compared to 
the values obtained for our SHARDS emitters. Our es¬ 
timates of the EWrf ([OH]) are consistent with those de¬ 
rived by other authors at similar redshifts (compare them 
with those from Kornei et al. 2012 and Teplitz et al. 
2003 in the figure), meaning that our technique based on 
medium-band data does not include any systematic bias 
on these measurements and provides robust EW values. 

The EWrf ([OH]) evolution curve appears to rise as 
—(l-bz)® up to z—1, and then flattens or starts to de¬ 
crease at higher redshifts. We interpret this as a real 
decrement presented by the EWrf ([OH]) evolution curve. 
In fact, the global behaviour of the evolution of the 
EWrf ([OH]) with redshift appears to follow that of the 
SFR density (SFRD) of the Universe (see e.g. Madau & 
Dickinson 2014, for a recent review), presenting a steep 
rise (proportional to —(l-|-z)“, with a—3.2), up to z—1 
and a flattening (or slow decline) at higher redshift (as 
suggested by the measurements at z—1.5 from Ly et al. 
(2012)). Indeed, the value of the slope a is quite sim¬ 
ilar to those derived form SFRD studies, where values 
between 3 and 4 are usually obtainded (e.g. Lilly et al. 
1996; Baldry & Glazebrook 2003; Perez-Gonzalez et al. 
2005; Tresse et al. 2007, 2002; Gucciati et al. 2012). This 
trend reflects the varying SF efficiency at different epochs 
in the life of the Universe, which is well traced by such a 
simple quantity as the EW of [OH] emitters. 

4. PHYSICAL PROPERTIES OF THE [OH] 
EMITTERS AT z-d.S^ AND z^l.23 

In this Section, we discuss and compare the SFR, stel¬ 
lar mass and (mass-weighted) age properties of the [OH] 
ELGs selected in the redshift bins identified by our two 
SHARDS filters. In the following Section we will com¬ 
pare the properties of [OH] emitters with the other sub¬ 
samples defined in Section 2.4, investigating the effect of 
using different SFG selection techniques. 

The stellar mass and the [OH] observed (not corrected 
for dust extinction) luminosity distributions, L([OH])obs, 
for the two samples of SHARDS [OH] emitters at redshift 
^ 0.84 and ^ 1.23 are shown in Figure 12 (left panels). 
The average values of the stellar mass in the two redshift 
bins are similar, passing from (log(M/Msun)) ~ 9.3 at 
z~0.84 to (log(M/Msun)) ~ 9.5 at Z'^1.23 (left panels). 
We also note that the average stellar mass is very simi¬ 
lar to that obtained for the parent samples for the two 
redshift intervals (overplotted as black histogram). 

In the right panels of Fig. 12 we show instead the dis¬ 
tribution of the observed and extinction-corrected [OH] 
luminosities (L([OH])obs and L([OH])cor), adopting the 
extinction values derived from our two-population SED- 
htting technique. The average optical extinction values, 
for the old and young populations, are Nom ^^0.9 and 


SHARDS ELGs at z-0.84 and z-1.23 


15 



Fig. 10.— Distributions of EWs, line fluxes and R-band magnitude for the [OII]-selected galaxies at 2:~0.84 (upper panels) and z~1.23 
(bottom panels). The vertical dashed lines represent the median values also given in each panel. The blue dot-dashed line in the left and 
central panels marks the minimum EW^f ([Oil]) and corresponding flux achievable assuming the typical uncertainty for the faintest galaxies 
in our sample. 
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Fig. 11.— Evolution of the median EWrf([OII]) with redshift. 
We represent the values derived in this study for 2 = 0.84 and z = 
1.23 and a compilation from the literature at different redshifts up 
to (see the legend). The dashed line represents the best fit in 

the form EWj-f ([OII])'^ (H-; 2 )“, with a = 3.2. The arrows indicate 
upper or lower limits respectively (see text for details). The datum 
at zr^O.Tb from Teplitz et al. (2003) should be considered as an 
upper limit since their detection limit does not allow to measure 
EWs lower than ~ 35A. The point from Kornei et al. (2012) should 
be considered as a (slightly) lower limit since their EW distribution 
is biased toward higher mass objects (<M>~1.5xlO^®Msun)- The 
data from Ly et al. (2012) show a lower value (^^ 50A) of the 
median EW measured for their two redshift bins at 2 :'^ 1.47 and 
2 :1.62 with respect to the trend presented in the lower redshift 
bins. These last points are not included in the fit, because they 
may represent a flattening of the EWr.f([OII]) evolution at 2 > 1. 



Fig. 12.— Left panels: distribution of stellar masses for the [Oil] 
emitters (hatched histogram) and the parent sample (black out¬ 
lined histogram). Vertical dashed lines indicate the median value 
for the [Oil] (red) and the comparison (black) sample respectively. 
Right panels: observed (hatched rose histogram) and extinction- 
corrected (hatched grey histogram) L([OII]) luminosities. Vertical 
dashed lines indicate the median value for the observed (red) and 
extinction-corrected (black) luminosity distributions. 


you '''1.5, respectively. We adopt the extinction from 


A, 

the young population to correct the [Oil] luminosities, 
since [Oil] emission is expected to be mostly produced 
in young star-forming regions. A Calzetti et al. (2000) at¬ 
tenuation law is assumed, with a factor ''.'2.3 accounting 
for the conversion from the stellar to the gas extinction. 
This conversion factor has been derived for a sample lo¬ 
cal star-forming galaxies, while recent studies have sug¬ 
gested that it might be evolving with redshift (see, e.,g., 
Wuyts et al. 2011; Kashino et al. 2013; Paimella et al. 
2015; Talia et al. 2015). Nonetheless, due to the large 
uncertainties in the derivation of this factor we preferred 
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Fig. 13. — Relations between (mass-weighted) stellar ages, stellar 
masses, and EWj-f ([Oil]) for the z = 0.84 (left panels) and = 1.23 
[Oil] emitters. These quantities were derived from two-populations 
SED fitting models that assumed a Chabrier IMF, exponentially 
declining SFHs, and Calzetti et al. (2000) internal dust reddening. 
See text for details on the SED modeling. 


to use the classical value for local galaxies. 

The median values for L([OII])obs are ^10^° ® erg 
and erg s-^ for the F687W17 and F823W17 

samples respectively, that shift to erg s“^ and 

erg s“^ for L([OII])cor- The median values for 
the observed and corrected [Oil] luminosity distributions 
increase about ~0.2-0.3 dex from the lower to the higher 
redshift bin. 

In Figure 13 we plot the M*-age and M*-EWrf([Oil]) 
diagrams for the two samples of [OII]-selected ELGs at 
z^O.84 and ~1.23. Despite the large scatter and intrinsic 
uncertainties present in the diagrams, we can draw some 
conclusion about the relation of [Oil] equivalent widths, 
stellar populations ages and galaxy stellar masses. More 
specifically, if we divide our samples in subsamples con¬ 
taining galaxies above or below the median value of the 
EWrf([OII]), we see that the median value of the mass is 
^4 ("^2) times larger for galaxies in the low-EWrf ([Oil]) 
sample with respect to galaxies belonging to the high- 
EWrf([OII]) sample at z^O.84 (2^1.23). 

Regarding the ages of the stellar populations we do not 
see any specific trend related to the EW of the [OII]- 
selected galaxies. Instead, we see that galaxies with 
masses below the median mass value have on average 
younger ages, by a factor of ^2 (~4), with respect to the 
more massive galaxies for the 2'^0.84 (2~1.23) sample, 
as might be expected within the downsizing scenario of 
galaxy evolution, where massive galaxies form the bulk 
of their mass at an early stage of evolution. These re¬ 
sults are in qualitative agreement e.g. with those from 
Ly et al. (2012), although we find larger average stel¬ 
lar ages perhaps due to our fitting method including two 
populations (thus being less affected by recent bursts of 
SF). 

Finally, in Figure 14, we plot the SFR-M* diagram for 
our [Oil] samples. In these plots we correct (assuming 
for the dust attenuation the prescription detailed below) 
the SFR derived from the observed L([OII])obs through 


the Kewley et al. (2004) relation: 

SFRobs([OII]) = 

(1.4 ± 0.4) X 10-41L([OII]) [M,,„ yr-i] ^ ’ 

and scale the SFR from Salpeter to Chabrier IMF by 
dividing by 1.7. 

We derive the extinction from the ratio of total SFR 
to the observed SFRobs([OII]), where the total SFR is a 
combination of the obscured and unobscured contribu¬ 
tions represented by the UV and IR SFR estimates (or 
predicted values from best SED-fits when dealing with 
IR non-detections). We adopt here the parametrisation 
introduced by Iglesias-Paramo et al. (2006, see also Bell 
2003; Hirashita et al. 2003; Iglesias-Paramo et al. 2004; 
Moncelsi et al. 2011): 

SFR = SFRobs(UV) + (1 - r/) X SFRobs(IR) (6) 

which includes a correction factor {1 — rj) accounting for 
the non-negligible contribution to dust heating coming 
from old stellar populations and that would bias the SFR 
estimate towards higher values. We assume rj = 0.32 and 
r] = 0.09 depending on the IR luminosity being lower 
than lO^^Lsun and higher than lO^^Lgun respectively (see 
Moncelsi et al. 2011). 

The observed SFRs were estimated for all galaxies in 
the parent sample based on either mid- and far-IR data 
from Spitzer and Herschel, or from UV luminosities. 
GOODS-N has been observed with the deepest MIPS 
data in the sky, with a Str limit of 30/iJy (Perez-Gonzalez 
et al. 2005), and also very deep PACS and SPIRE obser¬ 
vations, with 5(7 limits of 1.7 mJy and 9 mJy at 100/rm 
and 250//m, respectively. These limits (the deepest being 
MIPS 24/rm data) correspond roughly to 5-10 Mgun yr~^ 
for galaxies at 1.0<z<1.5 (Perez-Gonzalez et al. 2005; 
Elbaz et al. 2011). 

The total infrared luminosity, L(IR[8-1000]), has been 
estimated from SED-fitting of the MIR-to-FIR photom¬ 
etry using the Chary & Elbaz (2001), Dale & Helou 
(2002), and Rieke et al. (2009) templates. The infrared 
luminosities are then converted into SFRs using the Ken- 
nicutt (1998) relation : 

SFRobs(IR) = 2.7 X 10-44l(IR) [M^^n yr"^] (7) 

normalized to a Chabrier IMF. In order to estimate SFRs 
for galaxies not detected in the mid-or far-IR (including 
the quiescent galaxy sample), we used UV luminosities 
at rest-frame 280 nm. These luminosities were converted 
into SFRs by applying Kennicutt (1998) equation nor¬ 
malized to a Chabrier IMF: 

SFRobs(UV) = 0.95 X 10-2®L(UV) [M,,„ yr'i] (8) 

The UV-based SFRs were corrected for attenuation fol¬ 
lowing the recipe from Meurer et al. (1999), based on 
the UV/IR-/3 relation (IRX-/3). The UV slope (3 for each 
galaxy was measured from the SEDs, interpolating the 
best fitting models between 150 and 300 nm. 

We compare our results with the SER-M* relation from 
Whitaker et al. (2012), that studied the evolution of 
the SFR-M* relation up to redshift 2.5 for a sample of 
> 20000 galaxies from NEWEIRM Medium-Band Sur¬ 
vey, shown as the orange hatched region. Comparing the 
linear fit to our data (dashed black lines in Eigure 14) to 
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Fig. 14.— SFR derived from the [Oil] line versus stellar mass relation at redshift z~0.84 (left panel) and 2~1.23 (right panel). Observed 
SFR([OII]) (grey filled dots) and extinction-corrected SFR([OII]) (blue open dots) are shown in the plot. The black dashed line is the fit 
to our data, with the scatter in the relation shown as cyan shaded area. For comparison, we also overplot, as red diamonds with error 
bars, the median and (first and third) quartiles for the four equally populated mass bins including 90% of data (i.e. clipping extremes). 
The orange shaded area indicate the SFR-M* relation (and its scatter) from Whitaker et al. (2012) for the two redshift bins. The green 
continuum and dotted lines represent the SFR-M* at redshift z~ 1 and its 68% confidence interval, respectively, from Elbaz et al. 2007. 
In the right panel, we also depict as a grey dashed line the SFR-M* relation from the left panel. 


the SFR-M* relation at the same redshift from Whitaker 
et al. (2012) we observe a slightly larger slope in our case. 
This can be due to the different selections applied in the 
two cases, in fact a larger slope is also found by Whitaker 
et al. (2012) when selecting only blue SF galaxies, which 
might be more comparable with our [Oil] selection. 

Indeed, we should mention that there could be metal- 
licity effects, in the sense that for higher masses we ex¬ 
pect higher metallicities, and then we can have larger 
L([OII]) for the same SFR. In this sense, our SFR([OII]) 
could be overestimated at large masses. 

As check, we have also derived SFR([OII]) using the 
mass-dependent calibration for the dust-extinction from 
Garn & Best (2010). Using this relation, we would re¬ 
cover a better agreement with the results from Whitaker 
et al. (2012), in terms of slope of the massive end of the 
SFR-M* diagram. 

On the other hand, we compare our SFR-M* relation 
with the one presented by Flbaz et al. (2007) for redshift 
1 SF galaxies finding very good agreement between 
their logarithmic slope of ~ 0.9 and our best linear fit 
providing a logarithmic slope of ^ 0.98 at z ^ 0.84 and 
^ 0.86 at z ^ 1.23, and a scatter of ~0.35 and ^0.38, 
respectively, at these two redshifts. 

As noted in Whitaker et al. (2012) a flattening of 
the SFR-M* relation moving toward higher redshifts has 
been observed in several previous works (see e.g. Noeske 
et al. 2007; Karim et al. 2011) and confirmed using our 
[OII]-selected sample of SF galaxies using SHARDS data. 
Nonetheless, as also commented in Section 1, various ef¬ 
fects can contribute to bias these results, leading to a flat¬ 


tening of the SFR-M* relation, especially for the higher 
redshift samples, where the mass completeness may rep¬ 
resent a major issue. 

5. GOMPARISON OF SFG AND QUIESCENT 
GALAXY SELECTIONS 

In this Section, we compare the SF properties (as de¬ 
rived from SED fitting) of the [OH]-selected samples with 
the complementary samples introduced in Section 2.4 
(nominally, the mass-deleted, the IR-detected, and the 
quiescent UVJ-selected samples). This allows us to in¬ 
vestigate the effect of different sample selections on the 
determination of the main sequence for SFGs and to 
highlight the possible biases introduced by each selec¬ 
tion. 

As recently discussed in various works (see, e.g., 
Rodighiero et al. 2014; Speagle et al. 2014), the deter¬ 
mination of the SFR-M* relation using various sample 
definitions, related to different SFR indicators, can lead 
to different determination for the slope of this funda¬ 
mental observed relation, with a steeper relations typ¬ 
ically found when using UV-selected galaxies, with re¬ 
spect, e.g., to the IR-detected galaxy samples. This is 
also partially attributed to possible Malmquist bias in¬ 
troduced when using IR-selected galaxies, more easily 
detected at increasing masses. 

Nonetheless, as shown in Figure 15, this dependence 
of the slope on the sample selection is not clear in our 
samples at redshift z^l. In fact, we find very similar log¬ 
arithmic slopes (in all cases ~1) for the SFR-M* relation 
of the general parent sample and the different subsam- 
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Fig. 15.— Star formation rate versus stellar mass diagrams (as derived from SED fitting) for the different sample selections indicated in the 
legend. The dashed lines, identified by the color corresponding to the symbols in the legend, represent the linear fits to the SFR-M* relation 
for each subsample. The dark-grey dashed and dotted lines give the the SFR-M* for galaxies with a constant log(sSFR [yr~^])=—9.7 or 
—9.5 respectively. See the text for details. 


pies selected on the basis of their IR and [Oil] emission 
(represented by the black, magenta and blue dashed lines 
in Figure 15). For completeness, we also plot the best- 
fit linear relation for the subsample of quiescent (UVJ- 
selected) galaxies, which relies well below the main se¬ 
quence for SF galaxies. 

The total and normalized distributions of M*, SFRs 
and sSFRs, for each subsample are given in Figure 16. 
We observe that, at z=0.84, the [Oil] emitters have a 
very similar stellar mass distribution as the parent sam¬ 
ple, maybe with a deficit at low masses (as can be seen 
from the normalized histograms in Figure 16). How¬ 
ever, many galaxies (aproximately 50%) below log(SFR 
[Msun yi'~^])=0.8, i.e., 5-6 Mgun are not [Oil] emit¬ 
ters or are not detected by our survey. Since the ex¬ 
pected SFR detection limits, derived from the Scr flux 
detection limit for our survey in each filter (blue dot¬ 
ted lines in the central panels of Figure 10), is log(SFR 
[Msun yr-i])-(-0.2)-(-0.7), i.e., ~0.6-0.2 M^un yr’^, 
we conclude that this fraction of undetected low-SF 
galaxies are not due to a detection threshold, but have 
their [Oil] emission prevented/quenched by some physi¬ 
cal mechanism, e.g. low excitation due to the weak SF 
activity. In the case of IR emitters, it is clear that they 
concentrate in the high-mass end of the population, and 
probably the detection limits, of ~5 Mgun yr~^ (shown 
magenta dotted in Figure 16), avoids detecting low SF 
galaxies. A similar result is seen at 2=1.23. 

The distributions of sSFR for the different samples are 
given in the right panels of Figure 16. Interestingly, there 
is good agreement between the global sSFR distribution 
of the parent, [OII]-selected, and IR-selected samples of 
galaxies ). This similarity of the sSFR distributions, is 
mirrored by the presence of very similar main sequences 


for the various samples (see Figure 15) with a quite con¬ 
stant slope. Other recent studies have argued the pres¬ 
ence of a knee in SFR-M* relation, but we do not see it 
in our samples (see, e.g., Whitaker et al. 2014; Lee et al. 
2015; ?). As expected, the UVJ-selected sample of qui¬ 
escent galaxies exhibits a distribution of sSFR with very 
small values with respect to the other sample selections. 

Observing the normalised distributions, shown in the 
lower panels of each row of plots in Figure 16, we can con¬ 
clude that there is not (or only mild) mass-selection bias 
in the case of [Oil] emitters, while the IR-detected and 
UVJ-selected samples are systematically shifted toward 
higher masses. The peak in the SFR distribution is only 
slightly enhanced for the [Oil] sample, with respect to the 
parent sample, whereas the peaks for the IR-detected and 
UVJ-selected samples are clearly shifted towards higher 
and lower SFRs, respectively. Similar conclusions apply 
to both redshift ranges. 

As already discussed in Section 2.4, the assumed 
sSFR cut at log(sSFR [yr“^])=—9.7 (shown as dark-grey 
dashed lines in Figure 15) lead to a good overlap between 
the UVJ selection criteria and a pure sSFR selection for 
the lower redshift sample, while a slightly higher sSFR 
cut (log(sSFR [yr“^])=—9.5, shown as dark-grey dotted 
line in the right panel of Figure 15) would better reconcile 
the two definitions (i.e. UVJ-selected and sSFR-based) 
of quiescent galaxies for the higher redshift sample. This 
difference might be related to a small evolution of the 
SF properties between these two redshift intervals (cor¬ 
responding to ^^1.5 Gyr in time). 

Indeed we remark that a relatively large fraction of 
galaxies selected as quiescent using these two alternative 
definitions are detected as [Oil] or IR emitters. In fact 
we find ~20% and ^^25% of the lower redshift galaxy 
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Fig. 16.— Distribution of M*, SFR and sSFR for the z~ 0.84 (top panels) and z~ 1.23 (bottom panels) samples. Top panels in each 
plot represent the general number distributions, bottom panels give the histograms normalised to the peak, for each sample indicated in 
the legend (reported only in in the stellar mass plots). 


sample to have a [Oil] or IR detection, for the UVJ- and 
sSFR-selected samples, respectively. These fractions in¬ 
crease to —30% and —40%, for the higher redshift sam¬ 
ple. These results show, on one side, that the UVJ- 
selection seems to be more effective in the selection of 
bona fide quiescent galaxies with respect to a simple cut 
in the sSFR, on the other that moving to higher redshifts 
the contamination from SF galaxies start to be very rel¬ 
evant. This result is in agreement with recent results on 
the study of quiescent galaxies selected in the UVJ plane 
(see e.g. Belli et al. 2015). 

6. DUST ATTENUATION PROPERTIES 

Previous studies have found that the dust attenua¬ 
tion is a strong function of stellar mass (e.g. Reddy et 
al. 2006, 2010; Pannella et al. 2009; Wuyts et al. 2011; 


Whitaker et al. 2012), where the most massive galaxies 
are more highly obscured. 

We parametrize the attenuation as the ratio of the to¬ 
tal SFR (see Section 4 for the detailed definition) over 
the observed SFRobs for each SF indicator (UV1600, 
UV2800, [Oil] and IR respectively). We find a simi¬ 
lar trend using the SHARDS samples of [Ollj-selected 
galaxies in Figure 17: more massive galaxies have larger 
SFR/SFRobs ratios. We depict as solid symbols IR- 
detected galaxies for which the determination of the total 
SFR is more reliable, while arrows represent galaxies se¬ 
lected on the basis of their [011] emission but lacking IR- 
detection which may lead to an underestimation of the 
total SFR. Indeed, by selecting [OH] ELGs we are pos¬ 
sibly biasing the estimated average attenuation values 
toward lower values, since galaxies suffering strong dust 
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Fig. 17. — Top panels: ratio between the total SFR and the observed SFRobs for UV(1600), UV(2800), [Oil] and IR SFR indicators 
versus the stellar mass (left panel) and the total SFR (right panel) at redshift ~ 0.84 (F687W17). Arrows represent lower limits for the 
galaxies not detected in the IR. Dashed lines represent the median values computed for IR-detected galaxies, only. Bottom panels: the 
same for galaxies at redshift ~1.23 (F823W17). 


attenuation should have their [Oil] emission reduced or 
non-detectable. 

The IR alone is a good estimator of the total SFR for 
most galaxies, but there is fraction of them which have 
non-negligible amounts of dust-free star formation, and 
that seems to be more significant for lower masses. There 
is approximately a 20% (25%) of galaxies detected in 
the IR for which the contribution to the total SFR from 
unobscured star formation is larger than 30% at redshift 
^0.84 (^1.23), most of them being less massive than 
10^° Msun, and with total SFR lower than 20 Msun yr”'^. 

An analogous behaviour for the attenuation is found 
with respect to the total SFR (right panels in Figure 17), 
where the trend appears even more evident and the re¬ 
lation tighter exhibiting a smaller scatter. Galaxies sus¬ 
taining higher SFRs also exhibit larger dust attenuation 
(coherently with earlier results for SFGs, e.g. Dominguez 
Sanchez et al. 2012, 2014). This characteristic may 
be related to the fact that more actively SFGs possi¬ 
bly hold higher fractions of dust sustaining this high 
rate of SF or could be ascribed to a higher efficiency in 
heating the dust in these galaxies. The median value 
of the ratio SFR/SFRobs, where SFRobs is estimated 
from the UV(1600), UV(2800), [Oil] or IR observed (i.e. 
not extinction-corrected) SFR respectively, is found to 


smoothly increase as the wavelength of the SFR indica¬ 
tor decreases. 

Interestingly, the trend of the M*-SFR/SFRobs rela¬ 
tion in the two redshift bins considered does not seem 
to vary signihcantly suggesting no (or weak) evolution of 
the dust attenuation properties between these two epochs 
(separated by ^1.5 Gyr in time). Indeed, the attenua¬ 
tion of the line and continuum photons as represented by 
the [Oil] and UV(2800) SF indicators are quite similar, 
again suggesting that [Oil] can be considered a reliable 
SF indicator. 

The results presented in Figure 17 also suggest that 
galaxies of different masses are dominated by different 
mechanism of SF. The SFR is mostly dominated by the 
IR emission, when this is detectable, especially for more 
massive galaxies. This is indicated also by the tightening 
of the SFR/SFR/ij around 1 for high masses. Massive 
galaxies may present violent and explosive processes of 
SF (traced by high dust fractions), while less massive 
galaxies would undergo more relaxed episodes of SF. We 
can interpret these results as indicative of highly unsta¬ 
ble SF processes dominating in the massive galaxies, e.g. 
driven by galaxy interactions/mergers or violent relax¬ 
ation in highly unstable (clumpy) disks. Less massive 
galaxies would instead be characterized by a more diffuse 
SF. The investigation of the relation between morphol- 
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ogy/environment and SF properties is beyond the scope 
of this paper and is deferred to a dedicated work. 

7. SUMMARY AND CONCLUSIONS 

We have presented a detailed analysis of a complete 
sample of SFGs at 2—0.84 and 2—1.23. The sample 
is a comprehensive compilation of sources selected with 
3 different SFR tracers: the [Oil] emission-line, the 
UV emission, and the MIR/FIR emission. Concerning 
[Oil] emitters, we have selected them down to very faint 
magnitudes (—26.5 mag) benefitting from the ultra-deep 
spectro-photometric data from the Survey for High-z Ab¬ 
sorption Red and Dead Sources (SHARDS, see Perez- 
Gonzalez et al. 2013). The sample of [011] emitters has 
been compared with the general population of stellar 
mass-selected galaxies (including from quiescent to star- 
burst galaxies) at the mentioned redshifts. The SF ac¬ 
tivity and extinction properties of the whole population 
of galaxies at 2—0.84 and 2—1.23 have been character¬ 
ized with specihc SFR estimations, as well as UVJ colors 
and dust emission measurements based on Spitzer and 
Herschel data. 

The SHARDS data, combined with a customized and 
refined method for selection of ELCs using its medium- 
band (— 150 A) filters and for the determining the con¬ 
tinuum on the basis of SED-fitting, are suitable for ELC 
studies and have allowed us to obtain the results sum¬ 
marised below: 

1 . we have developed a technique that optimally ex¬ 
ploits the medium-band characteristics of the ultra¬ 
deep SHARDS survey to perform for the first time 
a systematic and robust detection of emission lines 
in low-, intermediate-, and high-redshift galaxies, 
inspired by the traditional procedures based on 
narrow-band imaging; 

2 . we have demonstrated that the depth and image 
quality of our survey allow us to recover virtu¬ 
ally all ELCs which have been already confirmed 
by the deepest spectroscopic surveys carried out 
in the COODS-N held. We have shown that 
we are able to extend these spectroscopic stud¬ 
ies of ELCs, typically limited to R—24-25 mag, 
to fainter magnitudes (i?—26-27 mag) and detect 
strong [OH] emitters with faint continuum emis¬ 
sion (i.e., large EWs), as well as Ly-a emitters up 
to 2—6 (Rodriguez Espinosa et al. 2014); 

3. by combining all the SHARDS hlters data and ht- 
ting them to stellar population synthesis models, 
we have shown that we can robustly measure equiv¬ 
alent widths and fluxes for these emission lines 
down to EWrf([OH])-15-20A and F([OH])-2- 
3x 10“^®ergs“^cm“^ respectively. These values 
are very similar to what can be done with deep 
spectroscopy. But, spectroscopic surveys typically 
pre-select their targets based on broad-band mag¬ 
nitudes (RI<24-25 mag) and have limited com¬ 
pleteness because they can only put a limited num¬ 
ber of slits in every object, unless huge amount of 
observing time is granted. By using our medium- 
band selection technique, we get all emitters with¬ 
out any strong selection effect and down to fainter 
broad-band magnitudes. 


4. we have focused this paper on two redshift ranges. 
One based on the selection of galaxies with emis¬ 
sion lines lying within the F687W17 filter (cen¬ 
tral wavelength: 687 nm), and another sample se¬ 
lected with the F823W17 filter. For [OH] emit¬ 
ters, these wavelengths correspond to 2—0.84 and 
2—1.23. The use of the SHARDS-based photomet¬ 
ric redshifts allowed to extend (a factor — 3 in num¬ 
bers) the samples selected based only on spectro¬ 
scopic redshift confirmation available from litera¬ 
ture. We can robustly measure [OH] fluxes and 
EWs even for galaxies that typically do not have 
sufficiently high-quality spectra for these kind of 
measurements (including galaxies usually exhibit¬ 
ing [OH] emission but lacking continuum detec¬ 
tion) ; 

5. the full SHARDS dataset, spanning the optical 
range from 500 to 950 nm, has allowed a signifi¬ 
cant improvement (about one order of magnitude, 
reaching azphot ^ 0.5% at redshift 2 — 1 ) on the 
photometric redshift accuracy (see Perez-Gonzalez 
et al. 2013; Ferreras et al. 2014, and G. Barro et 
al. 2015, in preparation) with respect to the use 
of previously available broad-band data. The ac¬ 
curacy in the photo-z determination is mirrored in 
the high success rate (>90-95%) and low contam¬ 
ination levels (<5-10%) of the selected samples of 
[OH] emitters; 

6 . we have analysed the EWrf([OH]) distribution for 
the two samples of [ 011 ] emitters, finding robust ev¬ 
idences on the evolution of the average EWrf ([Oil]) 
with cosmic time at least up to redshift — 1 , fol¬ 
lowed by a possible flattening, in good agreement 
with what is found studying the SFRD evolution 
of the Universe. Thus, globally speaking, the [OH] 
emission is a good tracer of the SFR evolution, de¬ 
spite the complicated dependence on intrinsic phys¬ 
ical parameters such as metallicity and tempera¬ 
ture, among others; 

7. we have found that galaxies with low-EWrf([Oil]) 
(i.e. below the median EWrf ([Oil])) have on av¬ 
erage higher masses with respect to galaxies with 
high-EWrf ([Oil]) (by a factor —4 and —2 for the 
samples at 2—0.84 and 2—1.23, respectively). We 
do not see a similar trend with respect to the age 
of the stellar populations. Instead, we observe 
that galaxies with masses below the median mass 
value have on average younger ages, by a factor of 
—2 (—4), with respect to the more massive galax¬ 
ies for the 2—0.84 (2—1.23) sample. As might 
be expected from the downsizing scenario, massive 
galaxies seem to have formed the bulk of their stel¬ 
lar mass at an earlier epoch; 

8 . we derive the SFR([OH]) and compare the SFR- 
M* relation for our samples with previous results 
from literature. We find a steeper slope {a —0.9) of 
the relation for both redshift intervals with respect 
to Whitaker et al. (2012), but this is compatible 
with their results when selecting blue galaxies (or in 
the case where we adopt ad mass-dependent extinc¬ 
tion calibration). Additionally, our data suggest a 
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trend of the slope with redshift, with higher red- 
shift galaxies showing shallower slopes (despite the 
various caveats related to possible observational bi¬ 
ases). Indeed the slope of the SFR-M* relation for 
our [Oil] emitters is in very good agreement with 
the main sequence of SF galaxies defined in 
Elbaz et al. (2007). Also the scatter in the relation 
(~0.35-0.38 dex) is in agreement with previous de¬ 
terminations; 

9. we have compared the distributions of physical 
properties and the SFR-M* relation for the [OII]- 
and IR-detected samples with those of the general 
population of galaxies at the same redshift (par¬ 
ent sample), finding no evident bias introduced on 
the SFR-M* relation by any of the selections of 
SFGs. Quiescent galaxies (identified on the basis 
of the UVJ diagram or a sSFR cut), populate a 
separated locus in the SFR-M* diagram and follow 
an independent relation, although there is a large 
scatter in the data points and some outliers with 
higher SFRs; 

10. we have discussed the possible biases affecting the 
different samples, comparing the total and nor¬ 
malised distributions for the M*, SFRs and sSFRs, 
for the subsamples of [OII]-selected, IR-detected 
and UVJ-selected quiescent galaxies, in the two se¬ 
lected redshift bins. We find only a very mild trend 
in selecting more massive galaxies based on the de¬ 
tection of the [Oil] emission line. A stronger bias 
is found relatively to the M* and SFR distributions 
of IR-detected galaxies, and (as expected) for the 
quiescent sample. We ascribe the larger bias for 
the IR-detected samples of galaxies to their larger 
mass- and SFR-detection limits; 

11. we observe very similar sSFR distributions for the 
parent sample and the [Oil]- and IR-detected sam¬ 
ples of SF galaxies. This behaviour is in relation 
to the very similar SFR-M* relations found for the 
various samples. We dot not observe significant 
variations in the slope of the main sequence of 
SF galaxies at high masses (M* > lO^^Mgun), as 
claimed by other recent works; 

12. we observe that the sample of quiescent galaxies, 
independently from the definition based on UVJ 
diagram or a sSFR cut, can be contaminated by 
both [Oil] and IR-detected galaxies. We find that 
~20% and ~25% of the lower redshift galaxy sam¬ 
ple present [Oil] or IR detection, for the UVJ- and 
sSFR-selected samples respectively. These frac¬ 
tions increase to ~30% and ~40%, for the higher 
redshift sample. We deduce that the UVJ-selection 
seems to be more effective in the selection of bona 
fide quiescent galaxies with respect to a simple 
cut in the sSFR, on the other hand we note that 
moving to higher redshifts the contamination from 
SF galaxies start to be very relevant and that 
more careful (redshift dependent) selection criteria 
should be adopted; 

13. we observe that, the IR emission alone is a good 
estimator of the total SFR for most galaxies, but 


there is fraction of them which have non-negligible 
amounts of dust-free star formation, and that 
seems to be more significant for lower masses. 
Approximately 20% (25%) of IR-detected galaxies 
present a contribution to the total SFR from unob¬ 
scured star formation larger than 30% at redshift 
~0.84 (^1.23), most of them being less massive 
than 10^° Msun, and with total SFR lower than 
20 Msun yr“^ 

14. we have finally found that the dust attenuation 
(as indicated by the ratio of the total to the ob¬ 
served SFR for each SFR indicator) strongly cor¬ 
relates with the stellar mass and total SFR of the 
galaxies. This result might be related to a higher 
dust content for more massive galaxies or to an 
higher dust heating efficiency in these massive ob¬ 
jects. The dust attenuation correlates linearly with 
stellar mass, and the trend of the stellar mass-dust 
attenuation relation does not show significant dif¬ 
ferences between the two redshift ranges, suggest¬ 
ing no (or slow) evolution of dust attenuation prop¬ 
erties between these two epochs (^1.5 Gyr apart). 
These results also suggest that galaxies of differ¬ 
ent masses are dominated by different mechanism 
of SF, with massive galaxies exhibiting more vio¬ 
lent and explosive processes of SF (traced by the 
higher dust fractions), while less massive galaxies 
would undergo more moderate episodes of SF; 

We have shown the huge potential of SHARDS ultra¬ 
deep data in studying galaxy evolution using ELGs. 
We plan to extend the results presented in this work 
using the whole SHARDS dataset (including 25 filters 
and several detected lines) to constrain the evolution of 
the physical properties of SFGs in the redshift interval 
0.3 < z < 2, fully covering the epoch where strong evolu¬ 
tionary processes shaped the galaxy populations toward 
the present-day Universe. 
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